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Abstract 


We have analyzed the discovery circumstances of all Earth- 
approaching asteroids detected in the last twenty-four years. 

In particular, we have calculated topocentric angular velocities, 
opposition distance, geocentric and heliocentric distances, phase 
angle, and lunar phase at discovery in an effort to separate any 
selection effects between chance and purposeful (i.e., as the 
result of a systematic search) discoveries. Another motivation 
was the possibility of discerning useful clues how to search more 
efficiently for such objects. There are 60 minor planets in our 
sample. Our principal result is that the discovery of Earth- 
approaching asteroids is dominated by serendipity. Therefore, 
searching for them at the current relatively bright limits at 
less than a very high rate seems pointless. 
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I. FORMULATION 


In this papei' any asteroid whose perihelion distance is 
less than that of Mars's is termed an Earth-approaching object. 
This includes the Aten, Apollo, and Amor classes and a few Mars- 
cross ing minor planets too. We have investigated the discovery 
circumstances of such objects with the intent of correcting the 
existing searches for systematic effects or exploiting such 
effects (if they exist) in order to enhance future search 
efforts. Because of the random first appearance on the celes- 
tial sphere of this usually fast-moving subset of the minor 
planets, any systematic advantage would be of great assistance. 
We looked for correlations amongst the discovery geocentric dis- 
tance, heliocentric distance, phase angle, topocentric angular 
velocity, angular distance from opposition, lunar phase at the 
time of discovery, and seasonality variations of discovery 
rates. 

In order to compute much of this information one needs an 
orbital element set. Orbital element sets for this subset of 
asteroids vary widely in quality. Factors influencing the 
accuracy of the element set include the total observational arc 
(days, months, or years herein), the distribution of astrometric 
data over this time span, the recovery of "pre-discovery" obser- 
vations (again days to years for some members of this group), 
and the quality of the individual observations themselves. For 
some objects several element sets are now available and 
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represent various stages in the differential correction process. 
It was impossible to choose element sets of uniform Quality. 

Given our desire to perform the essential research of this paper, 
some compromise had to be reached. The detailed references to 
the particular element sets we employed are given in Table I. 

The longest running Earth-approaching asteroid search is 
that of Helin and Shoemaker, 1979, It dates from 1973. (Our 
own search is only four years old; Taff, 1981.) We have used 
the last eleven years and the preceding eleven too as the base 
for our sample. This includes 60 minor planets. Our aim in this 
case was to try and balance, as much as possible, the pre-search 
era with the post-search era. In this way we hoped to minimize 
long term effects having to do with the weather, observatory 
existence, distribution, and equipment, the essential observing 
techniques, and even astronomical fashion. Bringing the study up 
to date, (e.g. 1983VA, and so on) would only increase the number 
of objects without good element sets. Backdating further, say 
to include the Palomar-Leiden survey would add 1960VA = 2061, 
1959EH = 1980RB1 = 2629, 1959LM, and three P-L objects (4789, 
6344, 6743) for which the discovery circumstances were not 
published. Again we have made a debatable choice. 

The original pair of discovery observations and orbital 
element sets have been taken directlv from the appropriate Minor 
Planet Circulars or International Astronomical Union Circulars, 


( 5 ) 


The catalogues of Marsden and Bardwell, 1982a, b, proved very 
helpful in this regard as did a complete list of Atsns, Apollos, 
and Amors kindly supplied to us by Dr. E, P, Helin. 

In the next section we discuss the calculation of our 
derived quantities and in the following one we analyze it in 
several wavs. Our conclusions are in the last Section. 
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II. DATA CONSTRUCTION 


a) Angular Velocities 

Ideally the minor planet left a long trail, on a well timed 
plate, so that both ends of the trail could be independently 
measured and reduced. Slightly less ideal are two different 
photographs from the same night. In a few cases days intervened 
between the initial pair of observations. In even fewer instan- 
ces there was only a single data point. Except in the last 
instance or when the textual remarks in an lAU Circular provided 
detailed information, we calculated Aa, A5, and At in the 
obvious way. The topocentric angular speed *(u and position angle 
P were then computed from 

(Aa/At)cos 6 ~ 01 sin P 
A6/At = 0 ) cos P 

(0 is given in degrees/day (to the nearest hundredth) and P is 
given in degrees (to the nearest one) in Table I [columns 4 and 
5; column 2 gives the original source of the data (MPC or lAUC) 
and column 1 gives the provisional designation]. Those values 
computed from the longer time spans or other less reliable means 
are indicated by a colon in Table I. 

b) Opposition Distance and Phase Angle 
We used the geocentric (solar system barycentric post-1980) 
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equatorial rectangular coordinates of the Sun (Earth) to compute 
the solar right ascension and declination at the discovery time. 
Then we coupled this linearly interpolated information from The 
American Ephemeris and Astronomical Almanac (The Astronomical 
Almanac) with the discovery location to compute the angular dis- 
tance from the opposition at the instant of discovery. This is 
listed in the tenth column of Table I. The discovery location 
came from the element sets whose MFC source is given in the 
third column of Table I. We assumed that there were no perturb- 
ing influences acting over the time spans involved and that the 
element sets were representative of an osculating heliocentric 
location and velocity. 

Phase angles were calculated from the above information 
according to their definition. They are listed in the ninth 
column of Table I. The heliocentric (r) and geocentric (R) dis- 
covery distances occupy columns seven and eight of the Table. 

c) Lunar Phase and Seasonality 

We anticipated that most discoveries would occur near New 
Moon rather than near Full Moon. We assigned Pull Moon a phase 
of 0, First Quarter a phase of 0.25, New Moon a phase of 0.5, 
and Third Quarter a phase of 0.75. Prom the discovery date, the 
lunar synodic period/ and old issues of The American Ephemeris 
and Nautical Almanac (or The Astronomical Almanac post-1980) the 
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phases in the sixth column of Table I were deduced. Another 
thing we examined was the seasonality variation — both monthly 
and in seasons. As almost all of this information is present in 
the provisional designation, we have not listed the actual dis- 
covery dates. 
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Table I, Data on Earth-approaching Asteroids 
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Table I. Data on Earth-approaching Asteroids (Cont 
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Table I, Data on Earth-approaching Asteroids (Cont'd. 
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Table !• Data on Barth-approachiny Asteroids {Cont*d. 
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Table I, Data on Earth— approaching P^teroids (Cont'd 
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III. DATA ANALYSIS 
a) Angulai: Velocities 

We examined the angular velocity in three fashions. First 
we created the (binned) frequency of angular speeds^ see Pig. 1. 
A bin size of 0?25/d was used (e.q., 0“/d ^ w 0J25/d and so 
forth) up to 3J25/d. There were 4 minor planets with higher 
angular speeds. The average value of w over the sampJ.e is 
2?05/d with a standard deviation about the mean (including 
Sheppard's corrections) of 3?38/d. The lack of low speed Earth- 
approaching minor planets needs no explanation because topocen- 
tric angular speed is the principal method of discrimination for 
these objects. We have no deep insight into the form of the 
curve. Next we looked at the distribution of the position 
angles (Pig. 2), As you can see it is almost uniform except for 
a slight northeast quadrant deficiency. The mean position angle 
is 196° and the standard deviation about the mean (including 
Sheppard's corrections for the 45“ size bins) is 95°. The 
northeast quadrant deficiency is especially evident when 
quadrant binning is used. The frequency distribution is now 8, 
17, 18, 16 around the compass. Furthermore, there is no corre- 
lation ui<*i.ween P and discovery date. Lastly we plotted, on 
polar coordinate graph paper, an angular velocity scatter 
diagram. The most striking feature about it is that there are 
no striking features. 
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b) Opposition Distance and Phase Angle 


Let be the heliocentiric location of the asteroid, ^ 
the heliocentric location of the Barth, and ^ be the geocen- 
tric (» topocentric) location of the asteroid. Then the phase 
angle i}> and the opposition distance $ given by 

-A ’-A -A*-E 

cos <t> « — ' cos $ = p V 

''a^A ‘^A‘E 


Since ^ and R/:^ is usually much smaller than 

x*E or r^, it follows that 

cos » cos <> + (RA/rE)-sin 

Hence, for most discoveries near opposition i|) < 4>. Although 
there is a strong opposition effect {« ) in minor planet 

phase functions and an additional linear decrease of the phase 
function in magnitudes with increasing phase angle (beyond 
e 6°), the histogram in Fig. 3 does not support the hypothesis 
that Earth-approaching asteroids are preferentially discovered 
near opposition. An equal area plot on the celestial sphere 
centered at opposition reinforces this conclusion. The phase 
angle at discovery frequency distribution is shown as the dotted 
histogram in Fig. 3. Apparently once an Earth-approaching minor 
planet becomes bright enough, and most are found with V < IS"*, 
it is detected by whomever photographed it and knew what to look 
for. 
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The geocentric and heliocentric distances in Table I show 
that Earth-approaching asteroids tend to be found '"hen they 
approach the Earth- We assume that this is a consequence of 
their increased brightness and angular speed. There are unusual 
cases though (e.g.» 1981ET3). The distributions are shown in 
Figs. 4 and 5. 


c) Lunar Phase and Seasonality 

We'll discuss seasonality first; the counts are given in 
Table II. There is a noticeable lack of non-uniformity due to a 
pronounced deficiency of discoveries during April, May, June, 
and July. We feel that the obvious Northern hemisphere bias of 
observatory distribution coupled with some elementary solar 
system geometry accounts for this (e.g., opposition is low in 
the sky at the summer solstice and the galactic plane high). 

Not surprisingly Spring is the lowest return season too: Winter 

17, Spring 8, Summer 16, and Fall 19, The big jump in the 
number of September discoveries (which props up the Summer dis- 
covery counts) is probably a result of the geometrical factors 
mentioned above as well as observer bias due to the Spring 
drought, the short Summer nights, a few months rest, renewed 
enthusiasm, good Pall weather, and so on. 

The distribution of lunar phases is shown in Fig. 6. As we 
anticipated there is a dearth of discoveries near Pull Moon and 
a maximum near New Moon, There is no bias towards the waning 
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Table II. Monthly Discovery Variation 
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Third Quarter Moon relative to the waxing First Quarter Moon 
which is significant. Such an effect might have been predicted 
on the grounds that a well planned search would commence during 
the waning Third Quarter phase. For then, should an interesting 
.asteroid be found, the maximum amount of dark time would be 
available for further observations. Alternatively, one's enthu- 
siasm might be higher at the beginning of a search period then 
at the end, it is easier to work in the evenings than it is in 
the mornings, and one's eagerness and alertness would be 
expected to drop immediately after a discovery because experi- 
ence has taught us of the rareness of any find. 

d) Discovery Circumstances 

What scientific research were the majority of the discover- 
ers conducting when they exposed the photographic plate on which 
the minor planet's trail appeared? Most of them were not 
searching for Earth-approaching asteroids but were engaged in 
some other type of astronomy. In addition, some of the aster- 
oids of interest credited to the Helin and Shoemaker, 1979 
search were found while they were pursuing other studies. As an 
example 1981VA was found during a search for high inclination 
asteroids (Helin, et al., 1982). As a second illustration 
1982DB was found while they were engaged in cometary studies, 
(JPL press relese #990 4/7/82 MBM; see also Aviation Week & 

Space Technology for May 10, 1982, pg, 51 or Sky and Telescope 
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1 


63 , 455, 1982.) These can clearly not be credited to their 

* 

systematic Earth-approaching asteroid search efforts. As a 
third example we have i981ET3, a "UCAS" discovery of R. Bus. 

In an effort to ascertain what fraction of the Earth-approaching 
asteroids discovered by and credited to the Helin and Shoemaker 
systematic search but should not be so attributed, we have 
critically examined their less well publicized discoveries. 
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IV. SEARCH STRATEGY AND PAYOFF 


What do we know about the pirobability of occuirrence, based 
on the 55 element sets irefei’enced in Table I, of an asteroid of 
this type? The means and standard deviations about the means 
are shown in Table III. Clearly these object can come from most 
CiUarters of the celestial sphere but retrograde motion is ruled 
out (so far). Their speeds of approach (or recession!) are 
highly variable, no season or lunar phase is preferred, and so 
forth. The only certainty we have apart from this randomness is 
that they will be especially bright near opposition. The optimal 
search strategy for this situation is intuitivel'^ obvious and 
rigorous calculations (Taff, 1984) verify that it can be 
improved upon by less than 10%. Part and parcel of an efficient 
use of telescope time, photographic plates, human enthusiasm and 
fortitude, and so on is to search when the other constraints on 
the problem do not interfere. Thus, one should search near each 
New Moon and not in the Summer. June 23 and September 23 are 
equally unfavorable times of year except for a weather differ- 
ential. Certainly Full Moons, July 15, and September 1 are to 
be avoided. 

When are these minor planets actually found? Sixteen were 
found during the Summer (effectively one-quarter of the total of 
60). Twenty-five were found outside of the optimum lunar phase 
range of (0.33, 0.67). While this is a lot closer to 1/3 of the 
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Table III, Oj.‘bital Element Statistics 


Element 

Mean 

Standard Deviation 

Semi-major axis 

1.90 A.U, 

0.55 A.U. 

Eccentricity 

0.489 

0.163 

Argument of perihelion 

168" 

101" 

Inclination 

17J7 

14?8 

Longitude of the Node 

184" 

96" 
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total tha’i it is to 2/3 of it, seven of these were Summer 
events 1 Thus, thirty-four of sixty asteroids were found at most 
improbable times. Because we already know of the accidental 
discovery of at least three others (e.g., 1981ET3, 1981VA, and 
1982DB), it becomes increasingly clear that careful planning and 
the husbanding of resources have little to do with the discovery 
of these objects. Serendipity rules the skies. 

We conclude that purposeful organized searches with current 
techniques (e.g., specifically with respect to search rates and 
limiting magnitude) have not had and will nob have a high payoff 
in terms of discoveries. Sky coverage, awareness of the quarry, 
and very hard work will continue to yield unpredictably large 
numbers of discoveries if additional objects remain to be found. 
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Figure Captions 


Pig, 1. Histogram o£ the number of Earth-approaching asteroid.? 

as a function oC topocentric angular speed at time of 
discovery. 


Fig. 2. 

Fig. 3, 
Fig. 4. 

Pig. 5. 

Pig. 6. 


Same format as Fig, 1 for position angle of topocentric 
angular velocity. 

Same format as Pig, 1 for opposition distance. 

Same format as Pig. 1 for geocentric distance at time 
of discovery. 


Same format as Pig. 1 for heliocentric distance at time 
of discovery. 


Same format as Pig, 1 for lunar phase (New Moon = 0.5) 
at time of discovery. 
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